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ABSTRACT 



The massive star forming region W3 was observed with the faint object infrared 
camera for the SOFIA telescope (FORCAST) as part of the Short Science program. 
The 6.4, 6.6, 7.7, 19.7, 24.2, 31.5 and 37.1 ftm bandpasses were used to observe the 
| emission of Polycyclic Aromatic Hydrocarbon (PAH) molecules, Very Small Grains and 

1 Big Grains. Optical depth and color temperature maps of W3A show that IRS2 has 

blown a bubble devoid of gas and dust of ~0.05 pc radius. It is embedded in a dusty 
Q H - shell of ionized gas that contributes 40% of the total 24 /im emission of W3A. This 

q . dust component is mostly heated by far ultraviolet, rather than trapped Lya photons. 

This shell is itself surrounded by a thin (~0.01 pc) photodissociation region where PAHs 
show intense emission. The infrared spectral energy distribution (SED) of three different 
zones located at 8, 20 and 25" from IRS2, show that the peak of the SED shifts towards 
longer wavelengths, when moving away from the star. Adopting the stellar radiation 
field for these three positions, DUSTEM model fits to these SEDs yield a dust-to- gas 



> 

mass ratio in the ionized gas similar to that in the diffuse ISM. However, the ratio of the 
IR-to-UV opacity of the dust in the ionized shell is increased by a factor ~3 compared 



to the diffuse ISM. 
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Subject headings: dust, extinction - ISM : general - ISM 



1. Introduction 

At early stages of their lives, massive stars are deeply embedded in their native clouds. These 
stars process their environment with energetic radiation and stellar winds, creating HII regions 
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and stellar wind driven bubbles ([Weaver et al.l 119770 , Overpressure of these regions drives their 
expansion into the molecular cloud and the HII region changes from hypercompact to ultracompact 
to compact. Eventually, the ionized gas breaks out, creating a so-called champagne flow that 
disperses the molecular cloud. The interaction between massive stars and their natal cloud controls 
the star formation process by compressing gas clumps through shocks - triggering subsequent 
star formation -, photo-ionizing and photoevaporating envelopes around low mass protostars in 
the cluster - creating prominent globules, fingers, and proplyds - and dispersing gas in ionized 
and neutral gas flows. This evolution of the region is accompanied by a processing of the cloud 
material. Specifically, dust c oagulates and this has major influence on the absorption properties 



(jQssenkopf Henninglll994l ; lOrmel et al 



2011 
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The radiation from the powering massive stars is absorbed by the surrounding gas and dust 
cloud. Only long wavelength (mid-IR to cm) emission from the ionized gas, dust grains, and 
molecules, and molecular emission can escape and be used to study the properties of the newly 
formed stars and the early evolution of the region. In the 3-20 um region, strong a nd broad emission 
features at 3.3, 6.2, 7.7, 8.6 and 11.3 fim domi nate the spectra (jPeeters et al.ll2002l ) due to Polycyclic 



Aromatic Hydrocarbons (PAHs) molecules (jTielens 



20081 ). These broad emission features are 



perched on a mid-IR continuum due to transiently heated very small grains (VSG) whose nature is 
not well understood ([Desert et al.lll990l ). At longer wavelengths (^30 /jm), emission by big grains 
(BG) - in radiative equilibrium with the radiation field - take over ([Desert et al.lll990l ). Each of 
these emission components has been "calibrated" as a s tar formation rate indicators on regions 
of star formation in nearby galaxies (ICalzetti et al.l 120071 ) but these empirical validations have no 
general theoretical underpinning. 

Here, we present mid-IR observations of the deeply embedded OB star cluster in the W3 
main region obtained by FORCAST on SOFIA in order to understand the interrelationship of the 
different mid-IR emission components, their carriers and their evolution in regions of massive star 
formation. W3 main is a star formation region located tow ards the Perseus arm at a distance o f 
1.95±0.04 kpc, accurately measured using maser kinematics (IXu et al.ll2006l : lHachisuka et al.ll2006l ). 
Infrared and radio observations have detected a number of bright sources in th e W3 main cloud 



Harris Wynn- Williams! (|1976l ) ; I Wynn- Williams. Becklin. Neugebauerl (119721 ) ; iTieftrunk et al 



(jl997l ). These objects span a range in evolutionary history with deeply embedded protostars, 
ultracompact compact HII r egions, all assoc iated with a newly forming cluster containing some 



15 OB stars (jOiha et al.l 12004 ; iBik et al.ll20110. A large, embed ded population of intermediate- to 



low- mass stars coexists with the OB stars (jMegeath et al.lll996l ). In this paper we focus on W3A , 
a shell-like com pact HII region, powered by IRS2, a 2-3 Myr old 06.5V star (jOiha et al.l 12004k 
Bik et al.l l201ll ) that provides a classic example of a wind driven bubble HII region in an early 
stage. 
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Observations and data 



W3 main was observed by SOFIA during the Short Science flight series on December 4 and 
8, 2010, at an altitude of ~43000 feet. We acquired images at 6.4, 6.6, and 7.7 microns in single 
channel mode and 19.7, 24.2, 31.4, and 37.1 microns using the dichroic beamsplitter in dual channel 
mode. The field of view of the instrument is the same for all the images, ~3.4'x3.2', and the pixel 
size is 0.768". The point spread function has a FW HM of 2.8'' at wave lengths up to 20 fim and is 
diffraction limited at the longer wavelengths filters (jHerter et al.ll2012l ). We used a chop throw of 
~5' and a chop frequency of ~5 Hz. The integration time for each chop pa i r was ^30 s. The data 
were processed and calibrated using the pipeline described in lHerter et al.l ( 2012 ). The processed 
images were aligned and averaged using an iterative sigma rejection algorithm to produce images 
with effective integration times of 150 s for 6.4 and 6.6 /im, 270 s for 7.7 //m, 240 s for 19.7 //m, 150 
s for 24.2 and 31.4 //m, and 330 s for 37 fim. The scale size and ori entation of the resulting images 
have been checked against Spitzer IRAC images (jRuch et al.l 120071 ) and agree well. The estimated 
3o" uncertainty in the calibration due to variations in flat field, water vapor burden, and altitude is 
approximately ±20%. A false color image of the whole W3 main region is shown in Figure [Hand 
a zoom in on W3A is presented in Figure [2^,. 

The FORCAST photometry was checked against the W3-IRS5 spectrum obtaine d by the Short 



Wave length Spectrometer (SWS) on board of the Infrared Space Observatory (ISO) (|Boogert et al 



20081 ). W3 IRS5 is a point source in our images and is isolated even within the large SWS beam. 
At all wavelengths, the tw o agree to within th e absolute calibration error of 20%. W3A was also 
observed with SWS/ISO (jPeeters et al.ll2002l ). We integrated the emission over the FORCAST 
filters in the relevant apertures of the SWS spectrum, and found an excellent agreement in flux 
(~20%) at all wavelengths except for the 6.4 /xm band where the FORCAST flux is a factor 1.9 
larger. Given that SWS calibration of extended sources is quite challenging, we deem this agreement 
excellent. Strong [SIII] and [Nelll] lines are present in the SWS spectrum, but the flux contribution 
of these lines to the FORCAST 19.7, 31.5 and 37.1 bandpasses is less than 6%. We show cross cuts 
towards the North in the different FORCAST images in Figure 

W3A was also observed wi th ISOCAM-CVF. We use the Highly Processed Data Products 
(HPDP) from the ISO archive feoulanger et alJhood ). Comparison shows that the FORCAST 
fluxes, F6.4 and F7 .7, are dominated by t he PAH 6.2 /iin and 7.7 fim features and the Fq.q traces the 
underlying plateau (jPeeters et al.1 120021 ) . By subtracting the Fq.q contribution from the Fq^ and 
F7.7 bandpasses, we derive the "pure" PAH emission. 
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Results 



3.1. Observational analysis of W3A 

Figure [I] shows a false color image of W3 main, combining three mid-infrared bands tracing 
the PAHs and warm dust. All the well known, bright infrared and radio sources are visible in the 
map, including the deeply embedded protostars associated with W3 IRS5, the ultracompact HII 
regions W3B, W3C, and W3E, and the compact HII regions, W3A, W3B, W3D, and W3H. In all, 
some 10 regions powered by bright O or early B stars are discernable in the mid-infrared. The 
W3 OB cluster of sources is embedded in low-level PAH emission in the deep 8 /im IRAC images. 
However, rather than UV photons leaking out of the confines of the (ultra) compact sources, this 
emission may be powered by distributed B stars in a more extended cluster, postulated on the basis 
of the extended [CII] 157 //m emission (jHowe et al.lll99ll ). 



In the radio, the compact HII region, W3A, shows a clear shell -like structure to the north 
(jWynn- Williams. Becklin. fe Neugebauerl Il972l : iTieftrunk et al.1 Il997l ). At mid-IR wavelengths, 
W3A shows a very similar structure (Fig. [2]) and the -F19.4, -F24.2 brightness distribution tracks 
the radio emission well (cf . Figure ED • W3A seem to be a cl assic example of a w ind-driven bubble 
expanding and sweeping up the surrounding cloud material (jWeaver et al.lll977l ). The radio emis- 
sion and the -F19.4 and -F24.2 warm dust emission originate from the swept up dense ionized gas. 
This shell is surrounded by a dense PDR (at ~20" from IRS 2) traced by the PAH emission (Figure 
[3|). The sharpness and the large relative variations of the PAH emission implies that the PDR is 
very thin (~1"; 0.01 pc) and unresolved. The shell-like structure of this region produces strong 
limb-brightening and the residual -F19.4, -F24.2 emission inside the cavity and the PAH emission seen 
within the ionized zone reflect this geometry. The inner cavity is presumably fill ed with hot gas and 
devoid of dust. However, there is no evidence for X-ray emission f rom the gas (Home r et al.l l2002) 
possibly due to the effect of mass loading from embedded clumps (IRedman et al.lll998T). We deem 
the a lternative model - "leakage" of the hot gas into the surroundings (IHarper-Clark fc Murray 
2009J) - less likely in view of the thin, dense shell enwrapping the ionized shell in the IR images. 



We used the F31.5 to F37.1 flux ratio to derive color dust temperatures (Toust) & n d optical depths 



at 37 //m (T37), assuming a modified black body for the dust (k^ ~ u" with (3 = 1.8; lAbergel et al 



201 ll). Total dust c olumn density maps have been derived from T37 using the extinction curve taken 
from IPraine! ( 2003 ) with Ry = 4 (Fig. [2]). Results are summarized in Table 1. For Ry = 3.1 (5), 



the column density decreases (increases) by a factor of 0.92 (1.1). Color temperatures and optical 
depths derived using different filter combinations show substantially the same global structure 
but the absolute values increase by 20% (^24.2-^37.1) and 30% (i ? i9.4-i ? 24.2)- We prefer the -F31.5- 
F37.1 color combination because this emission is expected to originate from the same dust population 
(i.e. BG; cf., Figure [2]). We note that T£, us t peaks at the ionized gas shell rather than at the IRS2 
position, which likely reflects an increase in the relative importance of emission by the cold-dust 
in the PDR as compared to the warm-dust in the ionized-shell for sightlines traversing the cavity 
(see section [3T3|) . 
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The 6.2 and 7.7 /im PAH emission and the 19-24 /mi emission observed in W3A are asso- 
ciated with the neutral PDR and ionized gas, respe ctively. It is well established that th e PAH 
emission in compact HII regions peaks in the PDR (jTielens et al.lll993l ; iBerne et al.ll2009l ). The 
Spitzer/MIPSGAL & GLIMPSE surveys have shown that many HII regions present bright MIPS 
24 /im e mission associated with the ionized gas an d are surrounded by a shell of bright 8 /xm PAH 



emission ([Anderson et al 



2011 



Carey et al.ll2009l ). With this in mind, we recognize that the star 



W3 IRS2b is also surrounded by a shell of PAH emission (Figure [2|) while bright 19 & 24 /im emis- 
sion is present in the ionized gas. We thus surmise that thi s star is still em bedded in its own HII 
region, separate from W3A. Perusal of JHK images of W3A (|Bik et al.ll201ll ) also suggests that the 
emission from IRS2b is separate from that of the main HII region. IRS2 has already "swallowed" 
the HII region associated with IRS2a and IRS2c and likely the region around IRS2b is next on the 
menu. We consider that the ridge of bright emission towards the south of IRS2a is still a remnant 
of the earlier, isolated HII regions phase of this star. 

From the FORCAST data, the region associated with IRS2b (08V) has a 6 — 40 /im luminosity 
of 1 x 10 5 in good agreement with the estimate based on near-IR photometry and spectroscopy 



(Bik et al. 



2011 



. The whole W3A region has a 6 — 40 /im luminosity of 4.1 x 10 5 L . Subtracting 
IRS 2b, the obse rved 6 — 40 /xm luminosity of W 3A is ~3 x 10 5 L . Based on the spectral type of 
Bik et al.l (|201lh . and using iMartins et al.1 hood ), we estimate that IRS2 (06.5V), IRS2a (08.5V), 
and IRS2c (B0.5V) emit 1.6 x 10 5 L , 6.6 x 10 4 L , and 4 x 10 4 L Q , respectively. The total 
luminosity of the stars associated with W3A is thus ~3 x 10 5 L . We note that the FORC AST- 
derived luminosity is just a lower limit to the total infrared luminosity, because far-IR photometry 
is lacking. Nevertheless, it seems that the luminosity of the known OB stars is sufficient to explain 
the observations and no major contribution from embedded YSOs is required, consistent with 
Ojha et al.l ( 2004I ). The 6 — 40 /im luminosity associated with the ionized gas shell is 1.1 x 10 5 L . 
Hence, we estimate that the (radial) UV (absorption) optical depth of the dust in the ionized gas 
is, t uv = - ln(L 6 _4o/-t'*) - 1- 



3.2. Dust Spectral Energy Distribution 



To study the dust characteristics in W3A, three zones were selected to extract Spectral Energy 
Distributions (SED) in 6"x6" apertures: At 8" north of IRS2 inside the HII region; at 20" north of 
IRS2 in the PDR, and in the molecular cloud at 25" north of IRS2 (Figure [2]). These SEDs show a 
pronounced shift of the peak emission towards longer wavelengths when moving fr om the HII region 



to th e PDR and on to the molecular cloud. The DUSTEM model (version 3.8; ICompiegne et al 



2011 



was used to fit these SEDs, adopting the grain size distribution for the diffuse high galactic 
latitude model. We adopt an 06V t he closest to t he spectral type of the central star and used the 
appropriate solar meta llicity model (jKurucz II 19931 ) scaled to match the number of ionizing photons 
(jTieftrunk et al.lll997l ). modified to our adopted distance of 1.95 kpc. We have evaluated these 
stellar fluxes at the projected distances of the zones. For the HII region position, we included the 
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FUV and EUV stellar emission, while for the PDR and molecular cloud zones only FUV is included. 
In addition, the molecular cloud radiation field has been attenuated by a radial optical depth of 
Tjjy = 1. Finally, IR fluxes for each of the DUSTEM model components (PAH, VSG and BG) 
are adjusted to fit the observations. Resulting fits are shown in Figure [2] and derived parameters 
are summarized in Table 1. The results illustrate that the 31.5 and 37.1 /im emission originates 
from BG, which also dominate the 24 jim emission from the ionized shell. The 24 fim emission 
in the PDR has an important contribution from VSGs. The good agreement between the model 
results and observation for the long wavelength emission (i.e. BG), demonstrates that the adopted 
radiation fields - the main free parameters - are reasonable. 



3.3. Lyman alpha dust heating 



The dust temperature p eaks on the ionized gas shell and this is very reminiscent of Lya heating. 
From the radio observations ( Tieftrunk et al.lll997l ). we derive the Lya luminosity of 3.2 x 10 4 L Q or 
only 30% of the observed 6 — 40 /mi luminosity associated with the ionized gas shell. In general, the 
contribution of trapped Lya photons to the dust heat ing (r^) relative to the stellar radiation (r*), 
is approximately given by r a /T* = (r/Rs) 2 x 1/tjjv ( Tielens et al,ll2005l ) with r the distance from 
the star (0.12 pc), Rs the Stromgren radius of the region (0.2 pc) and tjjv ~1 the radial optical 
depth of the dust (see above). In the middle of the ionized gas shell, Lya heating contributes ~25% 
of the dust heating. 



We used Cloudy version 10.00 (jFerland et al.l 1 199a ) to model the W3A To ust distribution 
and study the heating due to trapped Lya photons. Assuming an 06 star with a luminosity of 
2 x 1O 5 L0, and a shell density of 5000 cm -3 , the model shows that only for tjjv ^0.2 (a factor of 5 
less than observed) dust heating is dominated by Lya and the temperature distribution tracks the 
density structure. Hence, both the observed luminosity and the observed optical depth imply that 
stellar radiation dominates the temperature distribution. We attribute the rather constant dust 
temperature in the ionized shell to the morphology of the region. Specifically, the absence of a rise 
of the dust temperature towards the central star reflects the presence of an inner dust-free cavity. 
The geometry of a dust-free cavity, surrounded by a dusty ionized gas shell and PDR (Figure [3b), 
will tend to give a rather constant color temperature with projected distance from the star, because 
the observed SED is a superposition of all emission components (ionized shell &: PDR) along the 
line of sight. 



3.4. Dust-to-gas ratio 



The H- nuclei column density was calculated from the radio (ITieftrunk et al.lll997l ). assuming a 
shell-like geometry, to be Nh = 7.8 x 10 21 cm~ 2 . The dust column density derived with DUSTEM 
results then in a dust-to-gas mass ratio of 0.012, close to the diffuse ISM value. In contrast, the 
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Table 1. 



Parameters 


H II 


PDR 


Mol. Cloud 




(0.08 pc) 


(0.2 pc) 


(0.25 pc) 


(ergs s^ 1 cm" 2 ) 


247 


27 


7 


Observational results 


737.1 


0.06 


0.24 


0.55 


N Dust (10- 6 g cm- 2 ) 


211 


830 


1921 


N H (10 21 cm- 2 ) 


7.8 






T co \ OT (K) 


72 


64 


50 


Dust-to-gas mass ratio 


0.012 






DUSTEM results 


T Dust (100 nm) (K) 


78 


55 


44 


Lpah/Ur (10- 2 ) 


0.97 


3.74 




A v 


2.1 


32 


61 


R v 


4.2 


4.0 


4.3 
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radial UV optical depth derived from the observed 6 — 40 //m l u mino sity (r^y=l, section 3.1) 
translates into a dust-to-gas mass ratio of 0.035 using the iDraind (|2003l ) Ry = 4.0 model. While 
this difference seems small, we consider it significant. To phrase this differently, adopting dust 
properties appropriate for the diffuse ISM, the radial absorption optical depth by dust in the ionized 
gas would be T\jy ~4 and only ~ 2% of the number of i onizing photons emitted by the central star 
would be used to ionize the gas (jPetrosian et al.1 119720 . This would imply an enormous ionizing 
photon luminosity that is incompatible with the observed properties of the ionizing stars in this HII 
region and the observed IR lumin osity of the dust. This is a well-known issue, dating back to the 
earliest IR studies of HII regions (jTielens & de Jonidll979l ). and it seems that the dust properties 
in the ionized gas have to be different from those in the diffuse ISM. This may reflect either the 
selective destruction of PAHs and VSGs in the ionized gas -which are important contributors to the 
FUV dust extinction- or more general ly an increase in grain size during the preceding molecular 
cloud core phase driven by coagulation ( Ossenkopf fe Henningl[l994 ; Pagani et alJboiO : Ormel et al 



2011 



)• 
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Fig. 1. — Color image of W3 main region, c olors are: blue 7.7 /xm, green 19.7 iim, and red 37.1 
nm. Sources h a ve be en identified following I Wvnn- Williams. Becklin. &: Neugebauerl ()1972l ) and 
Tieftrunk et al.l (jl997l ) . The location of the four embedded OB stars in W3A are indicated by 
crosses and are IRS2c, IRS2, IRS2a, and IRS2b (left to right). The rectangular box indicates the 
direction of the cross cut for which the photometry is shown in Figure [3l 
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Fig. 2. — Left: Spectral Energy Distributions derived from the FORCAST observations (red 
crosses) and ISOCAM-CVF spectra (red line) in the three positions marked with squares in 2a. For 
a clearer comparison, the contributions from the ionized gas emission lines have been subtracted 
from the CVF spectra. Best fit DUSTEM models are shown as black lines. The green lines show 
the separate contributions from PAH molecules, Very Small Grains, and Big Grains (peaking left 
to right). Right: a) False color image of the W3A compact HII region. The three squares indicate 
three positions where spectral energy distributions have been obtained, and crosses mark the posi- 
tion of IRS2c, IRS2, IRS2a and IRS2b from left to right, b) Color temperature map, the peak in 
emission is located at the north of IRS2 (section 3.1). c) Optical depth map at 37.1 ^m. 
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Fig. 3. — a) Cross cuts along the rectangle in Figure [Has a function of projected distance from 
IRS2. All values have been normalized to the maximum along the cross cut. F^ e + F\™ e (black); 
^19.7) -^24.2 (red); ^31.5 + ^37.1 (blue), b) Derived colo r temperature and opti cal depth at 37^m (see 
text for details), along with the 2 cm emission from iTieftrunk et al.l (|1997l ). c) Schematic plot of 
W3A, the relative emission from the PDR increases as line of sights are further away from IRS2. 



